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We present a way for implementing an n-qubit controlled-rotation gate with three-level super- 
conducting qubit systems in cavity QED. The two logical states of a qubit are represented by the 
two lowest levels of each system while a higher-energy level is used for the gate implementation. 
The method operates essentially by preparing a W state conditioned on the states of the control 
c 2 ^ ■ qubits, creating a single photon in the cavity mode, and then performing an arbitrary rotation on 

04 ' the states of the target qubit with assistance of the cavity photon. It is interesting to note that 

the basic operational steps for implementing the proposed gate do not increase with the number n 
. of qubits, and the gate operation time decreases as the number of qubits increases. This proposal 

is quite general, which can be applied to various types of superconducting devices in a cavity or 
coupled to a resonator. 
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P^. I. INTRODUCTION 



I Multiqubit controlled gates are of importance in constructing quantum computation circuits and quantum informa- 
. tion processing. A multiqubit controlled gate can in principle be decomposed into the single-qubit and two-qubit gates 
^ ' and thus can be built based on these elementary gates. However, when using the conventional gate-decomposition 
protocols to construct a multiqubit controlled gate [1-3] , the procedure usually becomes complicated as the number of 
qubits increases. This is because single-qubit and two-qubit gates, required for constructing a multiqubit controlled 
-H , gate, heavily depends on the number of qubits. Therefore, finding a more efficient way to implement multiqubit 
^ ■ controlled gates becomes important. 

During the past few years, based on cavity QED technique, several theoretical proposals for implementing an n- 
qubit controUed-phase gate with ion traps, superconducting qubits coupled to a resonator or atoms trapped in a cavity 
have been presented [4-10]. These previous works opened a new way for the physical implementation of multiqubit 
controUed-phase (or controUed-NOT) gates, which play a significant role in quantum information processing, such as 
quantum algorithms [11,12] and quantum error-correction protocols [13]. On the other hand, experimental realization 
of a three-qubit controUed-phase gate in NMR quantum systems has been reported [14]. Moreover, a three-qubit 
quantum gate in trapped ions has been experimentally demonstrated recently [15]. 

The existing proposals in [4-10] for realizing an n-qubit controUed-phase gate can not be extended to realize an 
n-qubit controlled-rotation gate (denoted as controlled- i? gate below). Note that multiqubit controlled- i? gates are 
useful in quantum information processing. For instance, they can be used to construct quantum circuits for: (i) general 
$_i . multiqubit gates [16], (ii) preparation of arbitary pure quantum states of multiple qubits [17], (iii) transformation of 
quantum states of multiple qubits [18], and (iv) quantum error correction [19]. In addition, multiqubit controlled- i? 
gates can be applied to construct quantum circuits for implementation of quantum algorithms [20] and quantum cloning 
[21], and so on. In this work, we will focus on how to realize multiqubit controlled- i? gates with superconducting 
qubit systems. As is well known, superconducting devices have appeared to be among the most promising candidates 
for building quantum information processors recently, due to their design flexibility, large-scale intergration, and 
compatibility to conventional electronics. 

We note that if an n-qubit controlled- i? gate is constructed by using the conventional gate-decomposition protocols, 
2" — 3 two-qubit controlled gates would be needed (for n > 3) [1]. Thus, assuming that realizing any two-qubit 
controlled gate requires one-step operation only, at least 2" — 3 steps of operations are required, which increase with 
the number n of qubits exponentially. For instance, by using the conventional gate-decomposing protocols, 29 basic 
operational steps are required to implement a five-qubit controlled- i? gate, and 61 basic operational steps are required 
to realize a six-qubit controlled- i? gate. 

In the following, we will propose a way for implementing an n-qubit controlled- i? gate with three-level superconduct- 
ing qubit systems in cavity QED. The method operates essentially based on this idea: prepare a W state conditional 
to the states of the control qubits, create a single photon in the cavity mode, and then perform an arbitrary rotation 
on the states of the target qubit with assistance of the cavity photon. The W state used for the gate implementation 
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FIG. 1: Schematic circuit of an ra-qubit controlled- _R gate. If and only if the n ~ 1 control qubits on the filled circles (qubits 
1,2, and n — 1) are all in the state |1), a unitary rotation R{9) is performed on the two logical states |0) and |1) of the 
target qubit (qubit n), by |0) cos6'|0) -l-sineil) and |1) -sin 6*10) -t-cos^ll). 



is defined as follows 



y/n — 



= 5^P,|1)«("-^)|2), (1) 



where Pz is the symmetry permutation operator for qubit systems (l,2,...,n — 1), ll)®*-" |2) denotes the 

totally symmetric state in which one of qubit systems (1, 2, n — 1) is in the state |2) while the remaining n — 2 qubit 
systems are in the state |1) . For instance, the W state is (|112) + |121) + |211)) for n — 1 = 3. Note that W-class 

entangled states were originally proposed by Diir W et al. [22], which are useful in quantum information. 

As shown below, our method only needs 7 steps of operations, which is independent of the number n of qubits. 
Thus, when compared with the gate operations required by the conventional gate-decomposition protocols, the gate 
operations in this proposal are greatly simplified, especially when the number n of qubits is large. Furthermore, we 
note that the gate operation time for this proposal decreases as the number of qubits increases. The present proposal 
is quite general, which can be applied to various types of superconducting devices in a cavity or coupled to a resonator. 

This proposal requires adjustment of the level spacings of the qubit systems. For solid-state qubit systems such 
as superconducting devices, the level spacings can be rapidly adjtistcd (e.g., in 1 ~ 2 nanosecond timescale for 
superconducting qubits [23]), by varying the external parameters (e.g., the external magnetic flux for superconducting 
charge qubits, the flux bias or current bias in the case of superconducting phase qubits and flux qubits, see e.g. [24-28]). 
It should be mentioned that tuning the level spacings to have individual qubit systems coupled to or decoupled from the 
cavity mode was earlier proposed for the physical realization of quantum information processing with superconducting 
devices (e.g., see [5,29-31]). 

This paper is organized as follows. In Sec. II, we introduce the n-qubit controlled- i? gate studied in this work. 
In Sec. Ill, we discuss how to prepare the W state conditioned on the states of the control qubits in cavity QED. 
In Sec. IV, we present a way for implementing the n-qubit controlled- i? gate with three-level superconducting qubit 
systems in a cavity. In Sec. V, we give a brief discussion of the experimental feasibility for implementing a six-qubit 
controUed-Hadamard gate with superconducting devices coupled to a resonator. A concluding summary is provided 
in Sec. VI. 

II. A/^-QUBIT CONTROLLED-i? GATE 

For n qubits, there arc a total number of 2" computational basis states, which form a set of complete orthogonal 
bases in a 2"-dimensional Hilbert space of the n qubits. A quantum controlled- i? gate of n qubits considered in this 
paper is defined by the following transformation: 



iii2 ■ ■ ■ in-i) \in) ^ { ^'']' : ■ ^ Vr' T" n"-'=' " n" ' (2) 



l)i^(0)|^„), if nfc=^fc = i 
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FIG. 2: (Color online) Illustration of three-level superconducting qubit systems. In (a), the level spacing between the two 
upper levels is larger than that between the two lower levels. In (b) , vice versa. 



for all ii, ^2, • • •, in G {0, 1} . Here, the subscripts 1, 2, and n — 1 represent the n — 1 control qubits (1, 2, n — 1) 
while the subscript n represents the target qubit n; and \i1i2 ■ ■ ■ in-i) \in) is the n-qubit computational basis state. 
The operator R (0) is described by the following matrix 

^ /cos -sine\ 

^ ' \ smO cost) I ^ ' 

in a single-qubit computational subspace formed by the two logic states |0) ~ (l,0)"^and |1) = (0, 1)"^ of the target 
qubit n. It can be seen from Eq. (2) that if and only if the n — 1 control qubits (1, 2, n — 1) are all in the state 
|1) , a unitary rotation R{0) is performed on the two logic states |0) and |1) of the target qubit n. The definition of 
the n-qubit controlled- i? gate here is also shown in Fig. 1. 

III. PREPARATION OF THE W STATE CONDITIONED ON THE STATES OF THE CONTROLS 

In this section, we will discuss how to prepare the W state given in Eq. (1), when the n — 1 control qubits are 
initially in a computational basis state . For the gate purpose, the remaining 2"^^ — 1 computational states of 

the n — 1 control qubits need to be not affected during the W state preparation. In this section, we will also give a 
discussion on how this can be achieved. 

The superconducting qubit systems have the three levels shown in Fig. 2. Note that the three-level structure in 
Fig. 2(a) applies to superconducting charge-qubit or flux-qubit systems [24,25] and the one in Fig. 2(b) applies to 
phase-qubit systems [26,27]. In addition, the three-level structure in Fig. 2(b) is also available in atoms. In Fig. 2, 
the transition between the two lowest levels |0) and |1) is assumed to: (i) be forbidden due to the selection rules, 
(ii) very week due to the potential barrier between the two lowest levels, or (iii) highly detuned (decoupled) from the 
cavity mode during the gate operation. Throughout this paper, the two logic states of a qubit are represented by the 
two lowest levels |0) and |1) . 

To simplify our presentation, we will restrict our discussion to the three- level structure in Fig. 2(a). However, it 
should be mentioned that the results presented in this section and the method proposed in next section for the gate 
implementation are applicable to the quantum systems with the three-level structure depicted in Fig. 2(b). 

Consider n—l three- level superconducting qubit systems (1, 2, n — 1) in a single- mode cavity. The n—1 systems 
play the role of controls in realization of the n-qubit controlled- i? gate, discussed in next section. Assume that the 
cavity mode is coupled to the |1) O |2) transition but does not affect the level |0) (Fig. 3), which can be achieved by 
prior adjustment of the level spacings [24-28]. The Hamiltonian is given by (assuming h = 1) 

H = ujoSz + Wcfl+a + g + a5+) , (4) 

where 5, = i E"=i (|2>j (2| - (1|), 5+ = J^^^ll \2)^ (1| , = ^^=1 (2| , a+ and a are the photon creation and 
annihilation operators for the cavity mode, Wc is the cavity-mode frequency, wq is the |1) -o- |2) transition frequency, 
and g is the coupling constant between the cavity mode and the |1) o |2) transition. 



4 




|0) 

FIG. 3: (Color online) The transition between the two lowest levels is forbidden due to the selection rules, very weak due to the 
potential barrier between the two lowest levels, or highly detuned (decoupled) from the cavity mode. The |1) -f-*- |2) transition 
is non-resonantly coupled to the cavity mode with a detuning Ac, and non-resonantly coupled to the classical pulse with a 
detuning Ap. 



In the case when the detunning Ac — ujc — ojq ^ g^/rT+T with n being the mean photon number of the cavity 
mode, the Hamiltonian (4) can be rewritten as follows [32] 



H — uJoSz + uJcii^a — A 



^(|2)^.(2|-|l)^.(l|)a+a 



(5) 



where A = g^/Ac. The third term describes the photon- number dependent Stark shift while the last term describes the 
dipole coupling among the n—1 qubit systems. If the cavity mode is initially in the vacuum state |0), the Hamiltonian 
(5) reduces to 



(6) 



If the n—1 control qubits are initially in the computational basis state i.e., the Dicke state | J, — J) with 

J — [n — 1) /2, they evolve within the symmetric Dicke subspace spanned by {| J, — J) , | J, — J + 1) , | J, J)} . Here, 
the state | J, — J + fc) with k = 0, l,...,ri — 1 is a symmetric Dicke state with k systems being in the state |2) while 
n — k — 1 systems being in the state |1). The Dicke state | J, — J + k) is given by 



I J, - J + fc) = -jL= V |i)«("-'=-i) 12)^'= , (7) 

where Pz is the symmetry permutation operator for systems (1, 2, n— 1), ^ |i)®("~''~i) \2)'^^ denotes the totally 
symmetric state in which n — fc — 1 of systems (1, 2, are in the state |1) while the remaining fc systems are 
in the state |2) . 

One can check that the Hamiltonian i/o has the following properties 



Ho\J,-J + k)=ek\J.-J + k), (8) 

with 

=wo(-J + fc)_fc(2J-fc+l)A. (9) 

Eq. (8) demonstrates that the Dicke state |J, — J + fc) is the eigenstate of the Hamiltonian Hq with the eigenvalue 
Sk- The energy-level spacing between | J, — J + fc) and | J, — J -I- fc -I- 1) is e^+i — £k = ujq — 2 (J — k) X, depending on 
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FIG. 4: (Color online) (a) Non-identical level spacings for the energy levels. The level spacings between two adjacent levels 
become wider by 2A as the energy levels go up. (b) Illustration of the pluse (with frequency lo = ujq — 2JX) resonantly coupled 
to the transition between the two lowest Dicke states | J, — J) and | J, — J + 1) but detuned from the | J, — J + 1) -o- | J, — J + 2) 
transition with a detuning 2A. 



the excitation number of the state | J, — J + fc). It can be seen that the energy level spacings in the symmetric Dicke 
subspace are unequal [Fig. 4(a)]. For the detailed discussion, see Ref. [33]. 

To prepare the W state of Eq. (1), we now apply an external driving pulse (with frequency w) to the n~l systems 
(l,2,...,n— 1). Suppose that the pulse is coupled to the ]1) ]2) transition but far-off resonant with the transition 
between any other two levels of each system (Fig. 3). Thus, the interaction Hamiltonian between the pulse and the 
n — 1 systems is given by 



where is the Rabi frequency of the pulse. The Hamiltonian for the whole system is 

H = Hq + Hsp 



(10) 



(11) 



Performing the transformation U 



where 



^lutbi ^ obtain the engineered Hamiltonian in the symmetric Dicke subspace 



UHU+ - ujS^ 

n-l 

Y,h\J,-J + k) {J,-J + k\ 

k=0 
n-2 

+ J2^k{\J,-J + k + l} (J,-J+fcl +H.C.) 



(12) 



k=0 



rik = vt^/{2j -k) (fc + 1), 

5k = Wo (-J + fc) - fc (2J- fc + 1) A - cj (-J + fc) . 



(13) 



Assume that the applied pulse is resonant with the transition between the Dicke states ]J, — J) and \J,—J-\-l) 
[Fig. 4(b)]. Namely, the pulse frequency uj is set by w = — 2JA, i.e., Ap = wq — w = 2JA (Fig. 3). Discarding the 
constant energy — 2 J'^A we have 5o = = and 5k = k {k ~ \) X {k > 2). Hence the detuning of the pulse frequency 
with the transition frequency between the two energy levels [J, — J + 1) and \ J, — J + 2) is 2A [Fig. 4(b)]. Therefore, 
if we have (2J — 1) ^ 2A, which is guaranteed by setting 



fly/n - 1 < A, 



(14) 
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then the transition between the two Dicke states | J, — J + 1) and | J, — J + 2) is neghgible due to far-off resonance with 
the pulse. As a result, when the systems (l,2,...,n — 1) are initially in the Dicke state |J, — J) or \ J,—J-\- 1), the 
Dicke state | J, — J + 2) will not be excited by the pulse and therefore no transition from the state | J, — J + 2) to any 
one of the Dicke states {| J, — J + 3) , | J, — J + 4) , | J, J)} occurs. The Hamiltonian H' thus reduces to 

= rj\/2J (I J,-J+ 1) (J,-J| +H.C.). (15) 
It is straightforward to show from Eq. (15) that the states | J, — J) and | J, — J + 1) evolve as follows 

I J, - J) ^ cos(\/2JOt) I J, -J) - isin(\/2J0t) | J, -J + 1) , (16) 

I J, - J + 1) ^ cos(\/2jm) \J,-J+l)-i sin(\/2jm) | J, - J) . (17) 

Based on Eq. (1) and Eq. (7), it can be seen that the Dicke state | J, — J + 1) here is the W state defined in Eq. (1). 

From Eq. (16), it can be seen that when the control qubits (1,2, 1) are initially in the state | J, — J) (i.e., the 

computational state |ll..l)), the W state | J, — J + 1) is prepared through a transformation | J, —J) — >• —i | J, — J + 1) 
after a pulse duration t = 7r/(2-\/2Jf2). 

For the gate implementation below, we will need to transform the prepared W state | J, — J + 1) back to the state 
I J, — J) . Eq. (17) shows that this can be achieved by applying the same pulse to the systems (1, 2, n — 1) for a time 
t = 'k/{2s/2J^), via the transformation | J, — J + 1) — >• — 1 1 J, — J). 

The remaining (2"~^ — l) computational basis states \iii2---in-i) of the n— 1 control qubit systems can be classified 
into a set of Dicke states {|J — 1/2,— [J — 1/2))} with / = l,2,...,n— 1. For instance, the (n — l)-qubit computational 
basis states |11...110) and |11...100) can be written in term of the Dicke states | J — 1/2, —(J — 1/2)) (for I = 1) and 
I J — 1, —(J — 1)) (for I = 2), respectively. To see how the set of Dicke states {| J — 1/2,— [J — 1/2))} here not to be 
affected during the W state preparation, let us focus on a Dicke state \ J — 1/2, —(J — 1/2)) {I ^ 0), and discuss how 
to make this Dicke state unaffected by the pulse. 

The level spacing between the Dicke states | J — 1/2,— {J — 1/2)) and | J — 1/2,— {J — 1/2) + 1) is given by 

?! - eo = Wo - (2J - 0^- (18) 

Therefore, the detuning of the pulse frequency from the transition frequency between the two Dicke states 
I J — 1/2,— {J — 1/2)) and | J — 1/2,— {J — 1/2) + 1) would be w— (e'l— ?o) = —l^- The pulse Rabi frequency between the 
two Dicke states | J — 1/2,— {J — 1/2)) and | J — 1/2,— {J — 1/2) + 1) is fly/2J — I, which can be seen fr om the expres- 
sion of J7fe in Eq. (13) (for the present case, fc = and J is replaced by J — 1/2). If the large detuning J — I IX is 
met, the transition between the two Dicke states | J — 1/2, —(J — 1/2)) and | J — 1/2, —(J — 1/2) + 1) can be neglected 
due to far-off resonance with the pulse. Thus, the pulse does not excite the Dicke state | J — 1/2, —(J — 1/2) + 1) when 
the control qubits (1, 2, n — 1) are initially in the state | J — 1/2,— {J — 1/2)) . 

Note that for any I € {1,2, 1}, we have S}\/2J — I <C IX when Q^/n^^ <C A. Therefore, as long as the 

condition in Eq. (14) is satisfied, the Dicke state \ J — 1/2, —(J — 1/2)) with any given 1 — 1,2, or n — 1 will not be 
affected by the pulse, i.e., the rest (2"~^ — l) computational basis states of the control qubit systems (1, 2, n — 1) 
remain unchanged during the pulse. 

IV. IMPLEMENTATION OF AN A^-QUBIT CONTROLLED-i? GATE IN CAVITY QED 

Let us now consider n superconducting qubit systems (1, 2, n) placed in a cavity or coupled to a resonator. Each 
system has the three-level configuration. Initially, the transition between any two levels of each system is highly 
detuned (decoupled) from the cavity mode [Fig. 5(a, a')], which can be achieved via prior adjustment of the level 
spacings. In addition, assume that the cavity mode is initially in the vacuum state |0)^ . 

The procedure for implementing the n-qubit controlled- gate is listed as follows: 

Step (i): Leave the level structure of system n unchanged [Fig. 5(b')] while adjusting the level spacings of systems 
(1, 2, n—1) such that the |1) O |2) transition of each of systems (1, 2, n— 1) is non-resonantly coupled to the cavity 
mode, with a detuning Ac [Fig. 5(b)]. Then, apply a pulse to systems (1, 2, ...,n— 1) for a duration ti = -n / {2^/2JQ) 
[Fig. 5(b)]. As discussed in Sec. Ill, when the systems (l,2,...,n— 1) are initially in the computational basis state 
[11. ..1) , the W state \ J, — J + 1) is created after the pulse, via a transformation | J, — J) — >• —i \ J, —J + 1). Note that 
the cavity mode remains in the vacuum state during the operation of this step, as shown in Sec. III. 

Step (ii): Leave the level structure of system n unchanged [Fig. 5(c')] while adjusting the level spacings of systems 
(1, 2, n—1) such that the |1) o |2) transition of each of systems (1, 2, n — 1) is resonantly coupled to the cavity 
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FIG. 5; (Color online) The level structures of the systems (1,2, ...,n) during the gate preparation. Figures on the left side 
represent the level structures for systems (1, 2, ...n — 1), while figures on the right size represent the level structures of system 
n. Here, g is the non-resonantly-coupling constant between the cavity mode and the |1) ^ |2) transition, g' is the resonantly- 
coupling constant between the cavity mode and the |1) •f-> |2) transition, and g" is the resonantly-coupling constant between 
the cavity mode and the |0) -f-^- |2) transition. In addition, the transition between any two levels linked by a dashed line is 
highly detuned (decoupled) from the cavity mode and/or the pulse. 



mode for an interaction time t2 — (7i'/2) /{-/2Jg') [Fig. 5(c)]. The Hamiltonian describing this step is given by (in the 
interaction picture) 

Hi=g'aS+ +g'a+S-. (19) 

Here and below, g' is the resonantly-coupling constant between the cavity mode and the |1) O |2) transition. It can 
be found that under this Hamiltonian, the time evolution for the state | J, — J -f 1) ® |0)^ of the systems (l,2,...,n— 1) 
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and the cavity mode is described by 

I J, - J+ 1) |0)^ ^ cos{V2Jg't) \J,-J+1)<» |0), -ism{V2Jg't) \J, -J) (20) 

which shows that when systems (l,2,...,n— 1) are initially in the W state | J, — J + 1) , a single photon is created in the 
cavity mode after an interaction time t2 given above, through a transformation |J, — J + 1) ® |0)^ — > —i | J, — J) (8) . 
Note that the operation time t2 (oc 1/ y/n — 1) decreases as the number n — 1 of the qubit systems increases. 

Step (iii): Adjust the level spacings of systems (1, 2, n— 1) back to the original situation [Fig. 5(d)] such that the 
cavity mode does not couple to the systems (1, 2, n — 1). Meanwhile, adjust the level spacings of system n so that 
the |1) O |2) transition of this system is resonantly coupled to the cavity mode for an interaction time ts [Fig. 5(d')]. 
The Hamiltonian describing this step of operation is given by 

F/ = ;i(g'a+|l)(2|+h.c.). (21) 

The time evolution of the state [1),^ is described by 

|1)„ ^ cosig't) |1)„ -ism{g't) |2)„ |0), . (22) 

It can be seen from Eq. (22) that after an interaction time = 7r/(2,g'), the state [1)^^ changes to — i [2)^^ |0)^, . 
Note that the state |0)^ remains unchanged since the state |0)^ is not coupled to the cavity mode. Here and 
below, the subscript n represents the system n. 

Step (iv): Leave the level structure of systems (1, 2, n — 1) unchanged [Fig. 5(e)] while adjust the level spacings 
of system n so that the |0) O |2) transition of system n is resonantly coupled to the cavity mode for an interaction 
time t4 [Fig. 5(e')]. The Hamiltonian describing this step is 

Hi = h {g"a+ |0) (2| + h.c.) . (23) 

Here and below, g" is the resonantly-coupling constant between the cavity mode and the |0) O |2) transition. 
According to this Hamiltonian, one can easily find that after an interaction time ti = 0/g", the transformations 
|0)n |l)c ^ COS0 |0)„ - isin^ |2)„ |0)^ and |2)„ |0)^ ^ -isinO |0)„ + cosf? |2)„ |0)^ are obtained. 

The operations for the last three steps are the reverse operations of steps (i), (ii) and (iii) above, which are described 
below: 

Step (v): Leave the level structure of systems (1, 2, ...,n — 1) unchanged [Fig. 5(d)] while adjust the level spacings 
of system n such that the |1) O |2) transition of system n is resonant with the cavity mode for an interaction time 
[Fig. 5(d')]. The Hamiltonian describing this step is the one in Eq. (21). The time evolution of the state |2)^ |0)^ is 
described by 

|2)„ |0), ^ cosig't) |2)„ |0), -ism{g't) |1)„ |1), . (24) 

Thus, after an interaction time ^5 = 37r/(2,g'), the state |2)^|0)p becomes i|l)„|l)^. Note that the state |0)„|1)^ 
remains unchanged during this step of operation. 

Step (vi): Adjust the level spacings of system n such that the cavity mode does not couple to this system [Fig. 5(c')]. 
Meanwhile, adjust the level spacings of systems (1, 2, n — 1) such that the |1) o |2) transition of each of systems 
(l,2,...,n— 1) is resonant with the cavity mode for an interaction time te [Fig. 5(c)]. The Hamiltonian describing 
this step is the one in (19), from which one can easily find that the time evolution for the state | J, — J) (g) [1)^ of the 
systems (1,2, ...,n — 1) and the cavity mode is described by 

I J, - J) |1)^ cos{V2Jg't) I J, - J) O |1)^ -ism{V2Jg't) | J, - J + 1) (g) |0)^ . (25) 

It can be seen from Eq. (25) that the operation of this step results in the transformation | J, — J) (8) |1)^ — >■ —i \ J, — J + 1) 

(g |0), for te = Tr/{2V2Jg'). 

Step (vii): Leave the level structure of system n unchanged [Fig. 5(b')] while adjusting the level spacings of systems 
(1,2,..., n— 1) such that the |1) O |2) transition of each of systems (1,2,..., n— 1) is non-resonantly coupled to the cavity 
mode, with a detuning Ac [Fig. 5(b)]. Then, apply a pulse to systems (1, 2, n — 1) for a duration ^7 = t:/{2\/2JQ) 
[Fig. 5(b)]. As discussed in Sec. Ill, the operation of this step leads to the transformation | J, — J + 1) — >• — i | J, — J) . 

Note that after the last step of operation, we will need to leave the level structure of system n unchanged [Fig. 5(a')] 
while adjust the level spacings of systems (1,2, ...,n — 1) back to the original situation as shown in Fig. 5(a), such 
that the cavity mode does not couple to each system after the above manipulation. 
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The states of the whole system after each step of the above operations are summarized below: 

|11...1)|0)®|0), 
^*!^«_i|j,_j + l)|0)®|0), 

n'"^-|j,-j)|o)^|i), 

^*''-^"^-|J,-J)|0)^|l), 

Step^iv) _ (COS0 |0) - z Sine |2) |0)J 

Step(v) „ |_^^ _ (cosf? |0) + sin^^ |1)) 

1 1 J, - J + 1) (cosf? |0) + sine 1 1)) |0)^ 
step^vii) |ii__,i^(cog^|o) + sine|l))|0)^. (26) 



|11...1)|1)0|O), 
s*!^«_i|j,_j+l)|l)^|0)^ 

n^"'-|J,-J)|l)«|l), 
^"4"^|J,-J)|2)«|0), 
^'"^^^ I J, - J) (sin e |0> + 2 cose |2) |0) J 

Ste^(v) _ ^_ ^.^^ ^ cose 1 1)) |1), 

i I J, - J + 1) (- sin e |o) + cos e 1 1)) |0), 

step^vii) |n...i)(- sine |0)+ cose |1))|0)^. (27) 

where I J, — J) (i.e., |ll...l)) and | J, — J + 1) are the Dicke states of systems (1, 2, n — 1), while |0) , |1) , and |2) are 
the states of system n. 

On the other hand, it is noted that the following states of the whole system 

{\ili2...tn-l) |0)„ |0), , \hi2...ir^-l) |1)„ |0) J (28) 

(for all ii,i2, ■ ■ € {0, 1} and Hfe^i ^fc = 0) remain unchanged during the entire operation. This is because: (a) 

During the operation of step (i), the states {\iii2--.in-i)} of systems (1,2, ...,n — 1) were not affected by the applied 
pulse, as discussed in Sec. Ill; and (b) No photon was emitted to the cavity during the operation of step (ii), when 
systems (1,2, ...,n — 1) are in any one of the states {\iii2---in-i)} ■ Hence, it can be concluded from Eqs. (26) and 
(27) that the transformation (2), i.e, the n-qubit controlled- i? gate, was implemented with n systems (i.e., the n — 1 
control systems (1, 2, n — 1) and the target system n) after the above process. 

The systems not involved in each step of the operations above need to be decoupled from the cavity field and/or 
the pulse. This requirement can be achieved by adjusting the level spacings (e.g., doable for superconducting devices 
as discussed in the introduction) . 

The detunings Ac and Ap are set identical for each of systems (1,2, ...,n — 1) in steps (i) and (vii), and systems 
(1,2, ....n — 1) are brought to resonance with the cavity mode in steps (ii) and (vi). Therefore, the level spacings 
for systems (1,2, 1) can be synchronously adjusted via changing the common external parameters of the qubit 
systems during the entire operation. In addition, the cavity mode is virtually excited during the operation of steps 
(i) and (vii). Thus, decoherence caused by the cavity decay for these two steps is greatly reduced. 

From the description given above, it can be seen that the level |0) of each of qubit systems (1,2, ...,n — 1) is not 
affected during the entire operation, because the cavity mode was set to be highly detuned (decoupled) from the 
|0) o |1) transition and the |0) o |2) transition. Thus, the level spacing between the two levels |0) and |1) and 
the level spacing between the two levels |0) and |2) are both not required to be identical for each of qubit systems 
(1,2, ...,n — 1). However, as shown above, the level spacing between the two levels |1) and |2) needs to be identical 
for each of qubit systems (1,2, ...,n — 1). Note that for superconding qubit systems, it is difficult to have the level 
spacing between any two levels to be the same for each qubit system, but it is easy to have the level spacing between 
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certain two levels (i.e., the levels |1) and |2) for the present case) to be identical by adjusting device parameters or 
varying external parameters [34]. 

Finally, it should be mentioned that nonuniform of the device parameters for qubit systems (l,2....,n — 1) may 
cause the coupling strength g or g' (i.e., the coupling strength between the cavity mode and the |1) O |2) transition) 
not to be the same for each of qubit systems (l,2,...,n — 1). However, it is noted that for a superconducting qubit 
system, the qubit-cavity coupling strength is adjustable by varying the position of the qubit system in the cavity. 
Thus, by having the qubit systems (1,2, ...,n — 1) located at appropriate positions of the cavity, one can have the 
coupling strength g or g' to be identical for each of qubit systems (1, 2, n — 1). 

V. POSSIBLE EXPERIMENTAL REALIZATION 

In this section, we give a discussion on possible experimental implementations. For the method to work: 
(a) During the operation of step (i) or step (vii), the occupation probability pi of the Dicke state \J, — J + 2) due 
to the I J, — J + 1) -H- I J, — J + 2) transition induced by the pulse, and the occupation probability p2 of the Dicke state 
\J-l/2,~{J-l/2) + l) due to the \ J - l/2,-{J ~ 1/2)) ^ \J - l/2,-{J - 1/2) + I) transition induced by the pulse, 
given by [35] 



Pi ^ 

P2 



02 + A2/[2 (n-2)]' 

02 

02 + {ixf / [4 („ _ / _ 1)] 



02 

- 02 + AV [4 (n - 2)] ' ^^^^ 



need to be negligibly small in order to reduce the gate error. 

(b) According to the discussion in Sec. Ill, the following conditions need to be satisfied: 



,g<Ac, OVn - K .9 /Ac, Ap = (n - 1) .gVA^. (30) 

Note that these conditions can in principle be achieved because: (i) the Rabi frequency O can be adjusted by changing 
the intensity of the pulse, (ii) the detuning Ac can be adjusted by changing the jl) -H- |2) transition frequency wq, 
and (iii) the detuning Ap can be adjusted by varying the pulse frequency w. 
(c) The total operation time is given by 



T = 7r/(0V^rn:) + 7r/(5'V^r^) + 27r/5' + O/g" + 8r„, (31) 

which shows that for a given Q,^/n — 1, the r decreases as the number n of qubits increases. Here, Ta is the typical 
time required for adjusting the level spacings during each step. The r should be much shorter than the energy 
relaxation time and dephasing time 'y2p of the level )2) , such that decoherence, caused due to spontaneous decay 
and dephasing process of the qubit systems, is negligible during the operation. And, the t needs to be much shorter 
than the lifetime of the cavity photon, which is given by = Q/2itVc, such that the decay of the cavity photon 
can be neglected during the operation. Here, Q is the (loaded) quality factor of the cavity and is the cavity field 
frequency. To obtain these requirements, one can design the qubit systems to have sufficiently long energy relaxation 
time and dephasing time, such that r ^ l'2r^l2■p^ ''^'^'^ choose a high-Q cavity such that r <C . 

For the sake of definitiveness, let us consider the experimental possibility of realizing a six-qubit controUed- 
Hadamard gate (i.e., the controlled-i? gate for 9 = 7r/4), using six identical superconducting qubit systems coupled 
to a resonator [Fig. 6(a)]. Each qubit system could be a superconducting cliarge-qubit system [Fig. 6(b)], flux-qubit 
system [Fig. 6(c)], or flux-biased phase-qubit system [Fig. 6(d)]. As a rough estimate, assume g/2'iT ^ 220 MHz, which 
could be reached for a superconducting qubit system coupled to a one-dimensional standing-wave CPW (coplanar 
waveguide) transmission resonator [36]. With the choice of g',g" ~ 5, Ag ~ 10^, 0/27r ~ 1.1 MHz (i.e., A/0 ~ 20), 
and Ta ~ 1 ns, one has r ^ 0.2 /is, much shorter than min{7^"'^, 7^"'^} ^ 1 /is [26,37]. In addition, consider a resonator 
with frequency I'c ~ 3 GHz (e.g., Ref. [38]) and Q ~ 5 x 10^, we have ~ 2.7 /US, which is much longer than the 
operation time r here. Note that superconducting coplanar waveguide resonators with a quality factor Q > 10^ have 
been experimentally demonstrated [39]. 

For the choice of Ac ~ lOg and A/0 ~ 20 here, we have pi ~ 0.02 and p2 < 0.04, which can be further reduced by 
increasing the ratio Ac/g and A/0. How well this gate would work needs to be further investigated for each particular 
experimental set-up or implementation. However, we note that this requires a rather lengthy and complex analysis, 
which is beyond the scope of this theoretical work. 
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FIG. 6: (Color online) (a) Sketch of the setup for six superconducting qubit systems (red dots) and a (grey) standing-wave 
quasi-one-dimensional coplanar waveguide resonator. A is the wavelength of the resonator mode, and L is the length of the 
resonator. The two blue curved lines represent the standing wave magnetic field, which is in the ^-direction. Each qubit system 
(a red dot) could be a superconducting charge-qubit system shown in (b), flux-qubit system in (c), and flux- biased phase-qubit 
system in (d). The qubit systems are placed at locations where the magnetic fields are the same to obtain an identical coupling 
constant for each qubit system. The superconducting loop of each qubit system, which is a large square for (b) and (c) while 
a large circle for (d), is located in the plane of the resonator between the two lateral ground planes (i.e., the x-y plane). Ej is 
the Josephson junction energy (0.6 < a < 0.8) and Vg is the gate voltage. In addition, the external magnetic fiux $ applied to 
the superconducting loop for each qubit system is created by the magnetic field threading the superconducting loop. 



VI. CONCLUSION 



In summary, we have proposed a way for implementing an n-qubit controUed-rotation gate with three-level supercon- 
ducting qubit systems in cavity QED. This proposal requires seven steps of operation only, which is independent of the 
number n of qubits. In contrast, when the proposed gate is constructed by using the conventional gate-decomposing 
protocols, the basic operational steps increase with the number n of qubits exponentially. Thus, when the number n 
of qubits is large, the gate operation is significantly simplified by using the present proposal. In addition, as shown 
above, the gate operation time for this proposal decreases as the number of qubits increases. This proposal is quite 
general, which can be applied to various types of superconducting devices in a cavity or coupled to a resonator. 
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